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Abstract 


In this article, we explore the discontinuous Galerkin finite element method 
for two-parametric singularly perturbed convection-diffusion problems with 
a discontinuous source term. Due to the discontinuity in the source term, 
the problem typically shows a weak interior layer. Also, the presence of 
multiple perturbation parameters in the problem causes boundary layers 
on both sides of the boundary. In this work, we develop the nonsymmetric 
discontinuous Galerkin finite element method with interior penalties to 
handle the layer phenomenon. With the use of a typical Shishkin mesh, the 
domain is discretized, and a uniform error estimate is obtained. Numerical 


experiments are conducted to validate the theoretical conclusions. 
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1 Introduction 
In this article, we consider a two-parametric singularly perturbed problem 
with a discontinuous source term 


—eyu" + e9b(a)u + c(a)u=f, 2£E€Q=(0,1), 
u(0) = u(1) = 0, 


(1) 


where 0 < €1, €2 < 1. In addition, the coefficient functions, b(x) and c(z), 
are assumed to sufficiently smooth such that b(x) > bo > 0, c(a) > co > 0, 
and 

P= o(z) — 5e26"(z) 2 0 > 0, () 


where bo, co, and 7 are positive constants and the source function, f(x), 
has discontinuity of type J at some interior point say x = d. Moreover, we 
also assume that [f](d) < C, where [f](d) = f(d+) — f(d—). Problems of 
these types have several real-world applications, like transport phenomena 
in chemistry, biology, financial mathematics, and so on. Several numerical 
techniques are available in the literature on finite difference methods and 
finite element methods; for example, see [1, 10, 14] for the finite element 
method (FEM) and [3] for the finite difference method (FDM). Moreover, few 
articles [7, 4] are available that deal with two-parametric singularly perturbed 
boundary value problems, in which the authors have used various difference 
schemes to establish uniform convergence. In [11], authors have studied the 
superconvergence of non-symmetric interior penalty Galerkin FEM for two- 
parametric singularly perturbed boundary value problems. 

The FEM is one of the efficient numerical methods among all other nu- 
merical methods to compute numerical solutions of differential equations. 


Kadalbajoo and Yadaw [5] established a second-order convergence of FEM 
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in maximum norm on Shishkin mesh for two-parametric convection-diffusion 
problem. Brdar and Zarin [2] achieved uniform convergence of FEM in the 
energy norm on Bakhvalov mesh, while Zarin [12] has proved uniform con- 
vergence on exponentially graded mesh with FEM. In spite of these works, 
no one has considered the two-parametric convection-diffusion problem with 
discontinuous coefficients or source terms, which also predominately appears 


in the modeling of real-world problems. 


On the other way, the nonsymmetric interior penalty Galerkin (NIPG) 
method uses penalty parameters to enforce both the continuity of the solution 
and boundary conditions. Another interesting factor about the NIPG method 
is to handle the unstructured mesh, discontinuous coefficients, and local mesh 
refinements. Furthermore, an NIPG method can be made stable and coer- 
cive in the sense of convergence by taking some user-chosen penalty param- 
eters. Moreover, in other variants of the discontinuous Galerkin method like 
symmetric interior penalty Galerkin method and incomplete interior penalty 
Galerkin method, we have to choose penalty parameters, and it is a herculean 


task to establish its coercivity. 


Therefore, in this present work, we consider a two-parametric convection- 
diffusion problem with a discontinuous source term and establish a second- 
order convergence (superconvergence) of the NIPG method in the norm ||-||_, 


induced by the corresponding bilinear form. 


The article is arranged in the following way: In Section 2, we study the 
continuous problem and the bounds for different components of the exact 
solution. Section 3 is all about the bilinear form with respect to the NIPG 
method, its properties, and the existence and uniqueness of the weak solution. 
Construction of the Shishkin-type mesh, error analysis, and establishment 
of convergence results are introduced in Section 4. The numerical results 
and discussions are kept in Section 5 that which supports our theoretical 


conclusions. 


Note: Throughout this paper, we have taken C as a generic constant, 
which is free from perturbation parameters €1, €2, and the discretization 


parameter, N. 
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2 The continuous problem and bounds for the solutions 


This section is designed to discuss the solution of the continuous problem 
(1), their smooth and layer components, and their respective bounds. In 
general, the solution to the reduced problem (1) does not satisfy the boundary 
conditions. Therefore we have boundary and interior layers in the solution. 
To handle these boundary and interior layers, let pp and p; be two solutions 


of the characteristic equation: 
—e€,p? (x) + €gb(x) p(x) + c(x) = 0. 


Moreover, po(z) < 0 and p;(x) > 0 describe the boundary layers at « = 0 


and x = 1, respectively. Let 


A= oe po(z), Ar = Be pi(2), 


where 4; and A, are defined as 


—eb(x) — \/e3b(x)? + 4e1c(z) 


a a 26) 
a 2 2 
i) aaa é2b(x) + /€2)(z) + derc(x) 
2€[0,1] 24 


In that case, the problems are categorized on the basis of the relationship 


between the perturbation parameters: 


Case I. When €; < €2 = 1, then \; = O(e7') and A, = O(1). It is 
the case of singular perturbation problem with a only one perturbation pa- 


rameter. 


Case II. When €; < €2 < 1, then \; = O(ezey") and A, = O(e7'). 
Clearly one can see that ; < A;, which results a stronger boundary layer at 


x =1 than the boundary layer at « = 0. 


Case III. When <3 < €, <1, then ; and X,. both are of order Ole; ””), 


and problem behaves like reaction diffusion problems. 
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Further details about the category of the problems can be found in [11]. 
In this work, we discuss the a posteriori error estimation for the problem that 


falls under the category of case IJ. 
Theorem 1. Let p € (0,1) be arbitrary and let 
2m [8 l|z.o¢a) $B). 
Then, the following bounds for the derivatives of the exact solution holds: 


|u* (x) | < C(1 4 he Pe. Ape eS + Veet + eee. 


xrEQ, forO<k<m. 


Proof. The proof can be found in [6]. 


This motivates us to decompose the exact solution wu of (1) into smooth 


and layer components in the following way: 
u=S+ Ea he Fiz =F Ey “fF Ey 


such that the bounds satisfy the following relations: 


|S*(x)| < C, (3a) 
|EF,(@)| < CAFe?™*, (3b) 
|EF,(x) < Caper, (3c) 
|e oye e (3d) 
EF (ay) = Cake PA), (3e) 
forx€Qand0<k<m. 
3 The NIPG method 
Let Ty = J; = (#j-1,2;), 7 =1,2,--- ,N be a partition of the domain 


Q. Define the broken Sobolev space of order k, H fe (Q,Tn) : = {v : 
i7(O 2), 0], € H*(I;), for all I; € Ty}. Definitions of the Sobolev norm and 


seminorm are as follows: 
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N N 1/2 
2 
lellaay = (Dolley, ) toler = (Sob) 
j=l j=l 


We define the finite element space V™ that are related to the partition Ty 


1/2 


as follows: 


N= {v:veE L7(Q),v|, € P'(d;), for all I; € Ty}. 


I; 


Here, P1(J;) is the space of all polynomials of degree at most one on J;. It is 
to note that the elements in V™ are allowed to be discontinuous across the 


boundaries of each element J;. 


The corresponding weak formulation of (1) can be read as follows: Find 
un € VN such that 


B(un, vn) =Uun), for all uy € VN, (4) 


where 
B(u,v) = By(u,v) + Bo(u, v). 


Also, B,(u,v), Bo(u,v), and I(v) are defined for all u,v € H+(Q, Ty) as 


N 
Bu») => Jf atte Yoel ule] {v!(es)} ~ {u'(as) HC) 


(u,v) EI! e2d(x)ul + c(ax Wao Me v(x; ), 


Here, o;’s, j = 0,1,2,--- , N, are discontinuity-penalization parameters as- 


sociated with the grid points x;. Jumps and averages are defined by 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 347-366 


353 Discontinuous Galerkin approach for two-parametric convection-diffusion ... 


where v(t 


7) = lim, v(x) and v(x; ) = lim v(z), and at the boundary nodes 


=f = 


these are defined as 


[v(xo)] = v(#o),  {v(@o)} = v(20),  [v(aw)] = —v(ew),  {v(aw)} = v(2w). 


These notations can be seen in [8, 9, 15]. For any v € H?(Q, Ty), the norm 
is defined by 


N M N N 
1 
Jel = 21 oh So uye'(esg)? + Ihelldaay +Do (Foles) +09) ole? 
i=1 j=l j=l j=0 


Here, x;,; are the Gaussian points in the elements J; = (a;-1,x;), and w; > 0 
are the weights for M points Gaussian quadrature rule. The coercivity of 


bilinear form can be established by a simple calculation as in [8] as 
2 
Biv, v) = loll. 
Using coercivity, one can establish the result, 


llunlle $C llfllzaay + (5) 


where uy is the solution of weak formulation (4). More details can be found 
in [9]. Furthermore, the result (5) implies the uniqueness of the solution of 
(4). Rank-nullity theorem in a finite-dimensional space and uniqueness give 


the existence of the solution. 


Lemma 1. Let wu and uy be the solutions of (1) and (4), respectively. Then 


the bilinear form defined in (4) satisfies Galerkin orthogonality property 


B(u—uy,vy) =0, for all uy EV. (6) 


Proof. Detailed proof can be found in [11]. 


4 Error analysis on Shishkin-type mesh 


In this section, we introduce Shishkin-type mesh for the discretizing the do- 
main of definition. Then, we apply the NIPG method on the introduced 


mesh to avail the convergence result theoretically. Proceeding in this way, 
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the domain, 22 is divided into N subintervals as follows: Let Q11, Q2, Qs; 
Q,1, Qr2, and Q,,.3 be subdomains such that 


Qy4 = (0, do], Qi2 = (50, d — do]; Q7 3 = [d _ do, dl, Qr 1 = (d, d+ 51], 
Q,2 = [d+ 61,1—4)], OQ3 = [1 — 1, 1]. 


for some 69 and 6; defined by 


1l-d 
fo = min | $5 In} 61 =min{ >", = inv} 


Here, 6 is a user defined parameter. Therefore, the mesh length of each 


subdomain can be defined as 


WEN, for 111 UMN, 
ey for Qy 2, 
je (7) 
ee ln N, for Qy 1 U 0,3, 
cael I for Q,.2. 


The following bounds can be obtained for the above step sizes: 
CN-'InN);', for %1 UA 3, 
hps4 ON, for Q1,2 U O,.2, (8) 
CN-*InNd;z}, for Q,,1 UQ,3. 


Lemma 2. Let S/ be the piecewise Lagrange’s interpolation of smooth com- 


ponent S of the exact solution u of (1). Then, we have 
|s- si], son, 


Proof. Using the definition of jump for the smooth part of the solution, as 
S — S! is continuous across interelement boundaries. So, [(S — $7)(x;)] =0, 


it gives 


Now, we have to calculate the bounds for the two factors separately. We 


make the use of result on interpolation error estimate in [13] 
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N 


N M 
ed hi Y ((S — 81) (aig) i Cn LM ISli2c,) < cy 
i=1 20 j=l 


< C(N(N7'InNA;*)? + N74). (10) 


Here, we have used A>! < A; ' and bounds on the mesh lengths from (8). To 
establish the result in Lemma 2, we have to find out the essential bound for 


second term in (9) 


N N av; N 
Vile) las eO>5h if |S"(a)|\de < CS hi 
j=1 j=l j-1 


C(N(N7tInNA;')? +.N7?), (11) 


using (10) and (11) in (9). 


Lemma 3. Let Fy; and E,,; for 1 = 1,2 be the piecewise Lagrange’s inter- 
polations of left and right layer components FE; and E,, respectively. Then, 


we have the following estimation: 


||(Eui- a ae < CAN, i=1,2,s=0,1,2. 


||(Z,¢ — Et <CXN-*, i=1,2,s =0,1,2. 


Y lees = 


Proof. We will prove only the estimation for Ej, —E} ,, and other estimations 
for remaining part will be same. 


Making the use of the Cauchy—Schwarz inequality gives us 


(Bea — Bia) Mp1.) £ NBta) n= eyo + Ei) Me eyo.) 


S 211(F1,1) lay) - 


Here, we have used the stability property for special interpolant from Lemma 


3.2 in [15]. Now we have to bound only Il( E11)" le (aye, )" 


s s ALL s } SAT— 
(Fu) "Ile @a\anay S Or ae ONE SONS CALN 


Here we have used pd > 2, which is the required estimate, we had to prove. 


For estimations in the second part of the Lemma 2, we apply the same path 


of approach as above. So we have done the proof of lemma. 
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Theorem 2. Let E! denote the piecewise Lagrange’s interpolations of layer 
components of solution EF. Therefore, the interpolation error of the layer 


component satisfies 


|Z - E" ||, < C(N7* nN). 


Proof. To establish the bound in this theorem, we express errors in terms 
of left boundary layer, right boundary layer, and interior layer components; 
that is, 


|B - 2. < Bua — Ball +||B.2 — Eta). +l Bra — Frall-+|Br2— Erall.- 
(12) 

We derive the bounds for all the four terms in the right hand side of (12). 

Taking the first term from (12) and using coercivity of bilinear form from 


Section 3, we have 


N  M N 
2 2 
Ba — BIE =e oh (Ba BLY ead) +o ff ABB 
i=1  i=1 ga 
(13) 
Now, we estimate the first term from (13) over the fine mesh region 91, 


using classical result on interpolation error estimates 


N/S M P N/8 z; 2 
Ey S- hj x ((E11 = Ei 1) (#i3)) & Ce S- hi (/ ei lae 
i=1 j=l = 25-1 
N/8 ai XG 
< Ce y(f Evi Paz (/ a) 
j= Uj-1 Tj-1 
N/8 tee 
<Ce yf NeW 2PM da 
i= wget 
< C(N7'InN)*. (14) 
N/8 N/8 re 
Se Y(Eia — Ej,)’de < cS ny f |Ey,|?dx 
jai 75 j=l Uj-1 
N/8 


vj 
<CY A; | Ape “Ode 
j=l 8-1 


a 
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< C(N7'InN)*. (15) 


Using (14) and (15), we have established the first part of (12) for the finer 

part of the domain ();,;. In the same manner, we can establish the bound for 

Fia— Efal| , ||Bra — Bf lle and ||Er.2— Bal, over M2, Qa and Qa, 
||. i j ; 


respectively. Therefore, we can obtain 


2.2 - Efall.o,, SCN NY’, (16a) 
|| Ena — Fall etsy <O(N7'InN)?, (16b) 
||En2- Eralleo, $< CIN MN). (16c) 


Now, we estimate | #3 _ He | over .\ 07,1. In this domain, we will derive 
Tle 


the bound for the error in maximum norm. Using the result in Lemma 3 


N M N 
El b> h; > w; (Eri = Ei 1) (i,3)* < Ce, S- hyAZN~* < cn. 
i=N/8+1 j=l j=N/8+4+1 


Again, bound on the second term may be calculated using triangular inequal- 


ity and stability property as 


2 2 
|| = Fiallraoxa,s) = Fall e2caye,.) ae File @\any 


2 
<2 |Ei,1l 2 ¢a\e,1) 


<a} 


vj 
e€ 
j=N/8+41° 75-1 


Nee ONS, 


With the help of above two estimations a conclusion can be made 


2 


|| £1 — Era ones 


<CNn*. (17) 


Similar to (17), we can estimate the other remaining terms, and using (14)- 


(17), we conclude the bound in Theorem 2. 


Theorem 3. Let u’ be the piecewise Lagrange’s interpolation of the exact 
solution u of (1) on Gaussian points. Then interpolation error ¢ = u— ul 
satisfies 


IIc $ C(N* nN)’. 
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Proof. Using the solution decomposition and the estimations from Lemma 2 


and Theorem 2, we can conclude that 


Ile S$ |S-S']1.+ Be: — Bhall < O(N7! InN). 


k=l,r i=1 


Now, we conclude the bound for closeness error by means of Galerkin 


orthogonality and coercivity. Therefore, we need the following theorem. 


Theorem 4. Let u and uy be the exact and discrete solution of (1) and its 


weak formulation (4), respectively. For u/, the piecewise Lagrange’s interpo- 


T_ wy satisfies 


lation of u, let y =u 
IIxll. $< C(N™* nN)’. 
Proof. Coercivity and Galerkin orthogonality give us 


IIxll2 = Box, x) = —B(C,x). (18) 


From the definition, ¢ is continuous on interelement boundaries. That is, 
[¢(z;)] =0. So, Bi(¢,x) can be reduced to 


N a, N 
G0= af oxde- Valle ep- 9) 


Now, we estimate the first term in (19), 


N a; N vj 1/2, N 2; 1/2 
va | C'x'dx < (Saf (Car) (Saf (xr) 


j=l 


< C(N™ InN)? |Ix\le- 


Second term in (19) can be estimated using Cauchy—Schwarz inequality and 


the exact choice of discontinuity-penalization parameter a; = N, 


deca <(LSte oP) (Sainte?) © 


j=0 


< C(N* nN)? |Ixlle- 
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Using both these inequalities in (19), gives us 
Bi(¢,x) < C(N* InN)? Ix. - (20) 


Last but not the least is to estimate B2(¢, x), 


N ay N ny 
Df (cable! + e(a))xde = Sen [ole yenee 


We estimate the two terms of right side of above equation separately, 
N 
de 
j=l 
3 
<¢(SIKhme.o 
i=1 


< C(N™* InN)? |Ix\l.- 


| “ blr)Cy'da 


J 


3 
Xara + Dollz ean.) IX Urean0 | 
i=l 


Here, we have used 
Ixy S CA MN)? [Ix'Il2@,,) $ Cel Ar DN)? IIx. 


and results of Lemma 3. 
Finally, the last term can be bounded as 


Nan; 
S- / (cw) — eob(x))¢xda} < C [loll zeq IXllzaay < CCN? NY? Ilalle- 


j-1 


j=l 


Keeping in view of above two inequalities, we have 


Bo(¢,x) < C(N~* InN)? [Ix - (21) 


From (18), (20) and (21), we get the result of Theorem 4. 


Theorem 5. Let u and uy be the solutions of continuous problem (1) and 
the discrete problem (4), respectively. Then the discretization error satisfies 


the following estimates: 
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I|u — unl. < C(N7* nN)’. 


Proof. Applying triangular inequality on the estimation from Theorems 3 


and 4, we get the desired result. 


5 Numerical results and their assessments 


In this section, we consider two test problems to validate the theoretical 


estimation established in the previous section. 


Example 1. In this context, we have the first test problem: 


—e,u" + €2(3x? + 2)u’ + 7u = f(a), 
u(0) = u(1) = 0, 


where f(a) is defined by 


and d = 0.7. 


The above test problem is a singularly perturbed convection-diffusion 
problem with two parameters and discontinuity of type J in the source term 
at an interior point « = d of the domain. The solution of the above test 
problem exhibits boundary layers at both end points x = 0, 1. In addition, a 
layer at the point of discontinuity of source term x = d in the interior of the 
domain is recorded. The curve of the computed solution along with the exact 
solution is provided in Figure 1. The solutionprofile shows both the bound- 
ary layers at x = 0 and x = 1 and also the interior layer at x = 0.7. Here, 
0.7 is the point of discontinuity for f(a) in Example 1. The exact choice of 
discontinuity-penalization parameter used for this is 0; = N, for all 7, which 
are discussed in Section 4. 

Error and rate of convergence are examined for various values of dis- 
cretization parameter N and the singular perturbation parameter ¢,. Two 
parameters €, and €2 are related to each other by the relation ¢€, < ae For 


the sake of convenience in our computation, we take the general case €, = €3. 
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The exact solution to the above test problem is not known. Hence for the 
computational purpose, we determine errors for DG-norm by ||uy — u2n||pc- 
We apply the double mesh principle to deduce the rate of convergence and 


errors. Therefore, the following formula is applied 


RN =log, ( lun — van iipe ), 
7 ||uon = UAN IDG 


and the consequence is mentioned in Tables 1 and 2. From these tables, it 
is clear that the scheme is convergent with an order of convergence up to 
2. In the study of [11], we see that if we change our domain discretization 
from Shishkin mesh to Bakhvalov mesh and Bakhvalov mesh to exponential 
mesh, then the rate of convergence increases gradually. In this regard, we 
achieve comparatively better order of convergence, and numerically it will be 


achieved to 2. 


0.4 


0.3 4 


u(x) 


0.2 


0.1 


—— Computed Solution 


0047-77 Computed with 2N Grid Points 


0.0 0.2 0.4 0.6 0.8 1.0 
x 


Figure 1: Computed and exact solution for Example 1 for N = 2048 and ¢ = 1078. 


Table 1: ||jun — uen||pq errors for Example 1 


Number of intervals (JV) 
€2 64 128 256 512 1024 2048 
10-2 2.7829E-03 1.3574E-03 6.7059E-04 3.3331E-04 1.6617E-04 8.2964E-05 
10-%  3.7937E-03 1.8284E-03 8.9695E-04 4.4416E-04 2.2100E-04 1.1023E-04 
10-+  5.3400E-03 2.5120E-03  1.2140E-03 5.9633E-04 2.9548E-04 1.4707E-04 
10-° 7.7506E-03 3.5276E-03  1.6647E-03 8.0656E-04 3.9677E-04 1.9676E-04 
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Table 2: Convergence rates RN for Example 1 


Number of intervals (JV) 


€2 128 256 512 1024 2048 


10-? 0.8902 0.9406 0.9705 0.9855 0.9928 
10-3 1.1069 0.8868 0.9140 0.9545 0.9774 
10-4 1.7920 1.5795 1.0695 0.9100 0.9347 
i? 1.4870 1.8389 1.8514 1.5312 1.0518 


Example 2. The second test problem in this context is taken a 


—eyu" + €9(7a +5)u’ + 3u = f(z), 
u(0) = u(1) = 0, 


where f(a) is defined by 


and d = 0.7. 


The test problem given above in Example 2 is also regarded as two- 
parametric convection-diffusion problem with discontinuous source term. 
Due to the lack of an exact solution to the above problem, we use the double 
mesh principle to conclude the error and rate of convergence of the obtained 
numerical solutions. The rate of convergence is obtained by the following 


formula: 


PN = log, ( lun = u2n|lpg ) 
u2n — Uan|lpa 


The same choice of discontinuity penalization parameter as in the previous 
example is imposed here. The best possible relation between €; and €9 is 
€, = €3. Now, the error ||wy — u2n||pq and the rate of convergence P for 
the computed solution for the test problem in Example 2 with different values 


of discretization parameter N are given below in Tables 3 and 4, respectively. 
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The second-order convergence appears to be achieved up to the logarithmic 
factor from Table 4. 


Table 3: ||un — u2n||pq errors for Example 2 


Number of intervals (JV) 
€2 64 128 256 512 1024 2048 
10-2 9.8446E-03 3.0845E-03  9.9633E-04 3.3445E-04  1.1108E-04 —3.9958E-05 
10- 4.0731E-02 2.0132E-02 8.8730E-03 3.5715E-03 1.4107E-03  5.3149E-04 
10-4 4.1807E-02 2.0444E-02 8.8440E-03 3.5894E-03 1.3849E-03  5.3839E-04 
10-° 4.1841B-02 2.0351E-02 8.8344E-03 3.5877E-03 1.4237E-03 5.2625E-04 


Table 4: Convergence rates PN for Example 2 


Number of intervals (V) 


=) 128 256 512 1024 2048 


ig-* 1.6742 1.6303 1.5748 1.5901 1.4751 
10-2 1.8166 1.8819 1.6128 1.3401 1.4083 
iI ae 1.9320 1.8089 1.3009 1.3739 1.5630 
i" 1.9397 1.8239 1.6000 1.5333 1.5359 


5.1 Discussion 


Several articles are available in the literature that deal with the convergence 
and superconvergence of numerical solutions for problems with one and two 
small parameters. However, in all the articles, problems with continuous 
coefficients and source terms are discussed. In [11], authors have established 
the superconvergence result of the discontinuous Galerkin method for a two- 
parametric singular perturbation problem in which all the coefficients are 
assumed to be continuous. In that, authors have used three types of meshes 
usual Shishkin mesh, Bakhvalov mesh, and exponentially graded mesh for 


domain discretization and shown their rate of convergence of order two by 
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talking linear elements in its finite element space. The convergence results 
on both the Bakhvalov mesh and exponentially graded mesh have a sharper 
rate of convergence, that is O(N~?), than the rate of convergence obtained 
on the usual Shishkin mesh. In the present work, we attempt to obtain 
the convergence result for two-parametric problem in which the source term 
is discontinuous. The discontinuous Galerkin method is applied, and the 


second-order convergence has been obtained up to the logarithmic factor. 


6 Conclusion 


In this article, we proposed and discussed uniformly convergent discontinuous 
Galerkin FEM for two-parametric singularly perturbed problems with a dis- 
continuous source term. The problem usually exhibits a weak interior layer 
because of the discontinuous source term and the presence of the multiple 
perturbation parameters gives rise to boundary layers on both sides of the 
boundary. In order to address the layering phenomenon, we developed the 
non-symmetric discontinuous Galerkin FEM with interior penalties through 
this work. With the use of a typical Shishkin mesh, the domain is discretized, 
and a uniform error estimate was obtained. The numerical experiments were 


conducted to validate the theoretical conclusions. 
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